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SECTION  I 
INTRODUCTION 


Dynamic  plastic  response  of  metal  cylindrical  shells  subjected  to 
sharp  edge  blast  loads  is  mathematically  very  complex.  The  overall  plas- 
tic  deformation  process  is  complicated  by  the  buckling  phenomenon  associated 
with  the  conpressive  load.  IXie  to  the  highly  complicated  nature  of  the 
overall  response  many  approximate  analytical  solutions  coupled  with  exper- 
imental observations,  have  been  formulated.  The  general  results  of  typical 
approximate  solutions  (References  1,  2,  3)  are  expressions  for  a given  level 
of  damage  often  displayed  as  iso-damage  curves.  These  curves  are  generally 
a plot  of  impulse  versus  peak  blast  pressure  for  varying  values  of  the 
geometric  parameters  of  the  cylinder.  Considerable  experimental  work  using 
actual  blast  loads  by  Schuman  (References  4 and  5) , Lindberg  (Reference  3) 
and  presently  by  the  authors  has  not  completely  verified  the  approximate 
solutions.  / 

I so -damage  curves,  for  a given  damage  level,  may  be  drawn  using  the 
experimental  observations  and  the  approximate  methods,  but  a prediction  of 
failure  for  actual  material  separation  is  not  available.  Even  the  lar'e 
and  more  complex  computer  codes  (References  6,  7,  8)  are  limited  in  this 
area  due  to  the  lack  of  adequate  failure  criteria. 

The  main  objective  of  this  study  is  to  further  define  plastic  response 
of  cylindrical  shells  exposed  to  mild  transverse  blast  loads  which  produce 
material  rupture.  The  study  is  essentially  experimental  in  nature  with 
emphasis  on  failure  and  buckling  modes. 


SECTION  II 


EXPERINENTAL  TEST  ARRANGENENT 


Cylindrical  shells  of  a constant  radius  of  6.0  inches  of  varying 
lengths  and  thicknesses  were  tested  in  both  a fuel  air  explosion  (FAE)  and 
high  explosive  (HE)  environment. 

Testing  in  an  FAE  environment  was  accomplished  using  a gas  bag  tech- 
nique developed  previously  and  shown  schematically  in  Figure  1.  The  bag 
was  constructed  using  polyurethane  plastic  stretched  taut  on  a galvanized 
pipe  frame  and  held  together  with  3M  Paklon®  transparent  tape.  A 100-gram 
disc  of  green  Detasheet  and  detonator  were  placed  at  the  end  of  the  bag. 

The  bag  was  then  partially  filled  with  2 pounds  of  methyl  acetylene  propa- 
diene  (MAPP)  gas.  The  gas-air  mixture  was  then  mixed  using  a shaded  pole 
electric  motor. 

Detonation  of  the  Detasheet  produces  a Chapman- Jouget  wave  of  constant 
velocity  and  reflected  pressure  which  impinges  on  the  shell  at  the  opposite 
end  of  the  bag.  The  magnitudes  of  the  peak  pressure  and  impulse  were 
varied  by  adjusting  the  distance  between  the  end  of  the  bag  and  the  cylinder. 

Testing  of  several  cylinders  was  also  completed  in  an  HE  environment 
using  13.75-pound  50/50  Pentolite  spheres  hung  above  the  cylinder.  A 
schematic  of  this  test  is  shown  in  Figure  2.  The  charge  radius  (D/A)  was 
varied  to  give  peak  pressures  and  specific  impulse  similiar  to  that  of  the 
FAE  output.  Calibration  shots  of  the  spherical  devices  were  accomplished 
and  good  agreement  was  found  with  the  analytical  and  experimental  values 
given  by  Goodman  (Reference  9) . 
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SECTION  III 


BLAST  LOADING  ON  SHELLS 


’ 


To  analyze  the  response  of  cylindrical  shells  to  blast,  the  load  dis- 
tribution over  the  shell  surface  must  be  known.  Tor  the  purpose  of  this 
study  the  distribution  was  determined  experimentally,  and  is  reported 
here  for  possible  use  in  future  analytical  studies.  A thick  walled  steel 
cylinder  was  instrumented  at  its  mid  length  with  PCB  piezoelectric  trans- 
ducers. A series  of  tests  were  conducted,  using  the  bag  setup,  to  determine 
the  loading  from  the  fuel  air  explosive.  Figure  3 is  a schematic  of  the 
instrumented  shell,  and  Table  1 is  summary  of  the  results.  The  data  are 
presented  for  one-quarter  of  the  shell,  as  symmetry  is  maintained.  The 
data  in  Table  1 are  accurate  to  within  15  percent  or  better,  and  represents 
loading  estimates  on  shells  inside  and  outside  the  bag.  Data  for  a "IT' 
value  of  zero  represent  loading  on  the  shells  inside  the  gas  bag  from  the 
Chapman- Jouget  wave.  All  other  data,  for  D values,  of  3 to  6 feet,  repre- 
sent loads  from  an  air  shock  generated  from  the  bag  explosion.  The  peak 
pressure  was  curve  fitted  from  this  data  and  the  distribution  is  represented 
by  Equation  (1). 

P *(p  ♦ (P  - P )(cos  9)1-81  0°  < 6 < 90°  (1) 

m|_srs  j — ~ 

PR  * side-on  pressure  (9  3 90°) 

Pr  3 normal  reflected  pressure  (0  = 0°) 

The  peak  pressure  distribution  for  90  degrees  < 0 < 180  degrees,  or 
the  back  side  of  the  shell,  is  estimated  to  be  close  to  constant  at  the 
side-on  pressure  value.  The  complexity  of  actual  shock  interaction  and 
vortex  formations  on  the  back  side  is  beyond  the  scope  of  this  work.  A 
resonable  time  dependent  estimate  for  the  loading  may  be  given  by  Equation 
(2). 

P (0,t)  3 Pm[l  - t/T (eYJ  exp[-  a(0)t/x(9)]  ....0°  < 9 < 90°  (2) 

Note  that  u and  t are  functions  of  9,  and  must  be  obtained  from  the  data 
in  Table  1.  An  estimate  for  the  decay  constant  is  determined  by  first 
integrating  Equation  (2)  between  t equals  zero  and  t,  to  obtain  an  expres- 
sion for  the  cumulative  impulse  given  as  Equation  (3). 

1(9)  « Pm[T(9)/a(9)][l  * Va(9)  ♦ (l/a(9))(exp(-a(e)))J 


By  substituting  experimentally  determined  values  for  Pm  and  I from 

Table  1 into  Equation  (3),  a(6)  may  be  determined  for  each  value  of  D 

and  0.  If  the  shell  is  inside  the  fuel  air  mixture,  (D  ■ 0) , t is  reason- 
ably constant  as  seen  from  Table  1.  The  load  distribution  for  spherical 
charges  has  been  investigated  by  Zumwalt  and  Naik  (Reference  10),  and  this 
work  was  not  duplicated. 
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SECTION  IV 


EXPERIMENTAL  RESULTS 


A total  of  34  cylindrical  shells  were  subjected  to  blast  loads  from 
fuel  air  and  conventional  explosives  to  observe  plastic  buckling  and  fail- 
ure characteristics.  All  shells  were  made  of  6061-T6  aluminum.  Length 
was  varied  in  an  effort  to  determine  the  effects  of  these  parameters  on 
failure  and  buckling  modes.  The  magnitude  of  the  blast  load  was  controlled 
to  produce  material  rupture  or  near  rupture  whenever  possible. 

A suimary  of  the  test  results  appears  in  Table  2.  The  pressure  and 
impulse , P (0  = 0)  and  I (0=0),  are  those  reflected  values  as  seen  by  the 
leading  edge  of  the  shell.  The  normally  reflected  pressure  and  impulse 
values  for  the  explosive  sphere  tests  were  taken  from  Goodman's  report 
(Reference  9) . The  peak  pressures  are  accurate , but  the  cumulative  impulses 
may  be  slightly  in  error  due  to  the  relief  conditions  on  the  curved  shell. 

The  circumferential  mode  number,  n,  represents  the  whole  number  of 
buckled  waves  around  the  shell  circumference  assuming  the  cylinder  was 
fully  buckled.  The  number  n was  determined  by  counting  the  number  of  buckles 
and  dividing  by  the  fraction  of  the  circumference  that  was  buckled.  This 
nunber,  however,  is  not  truly  representative  of  the  deformation  process,  in 
that  only  a portion  of  the  shell  buckles  circumferentially.  The  percentage 
of  shell  deformed  was  relatively  constant  and  column  (11)  of  Table  2 shows 
that  only  25  to  35  percent  of  the  cylinder  circumference  was  plastically 
deformed. 

Before  examining  the  experimental  data  points  in  Table  2,  certain 
general  deformation  characteristics  should  be  mentioned.  For  all  cylinders 
tested,  a fundamental  mode  shape  was  observed  in  the  axial  direction;  i.e., 
m=l,  for  the  clamped  end  conditions.  Figure  4 shows  this  fundamental  shape 
for  a 0.071- inch  thick  shell.  Failure  in  the  cylinders  was  always  initiated 
at  the  clamped  boundary  and  approximately  at  the  0=0  point.  Rupture  of 
the  material  generally  propagated  from  this  point  in  both  directions  around 
the  cylinder  circumference;  and  the  resulting  failure  surface  resembled 
simple  tensile  failure  in  thin  sheet.  A typical  failure  is  shown  in  Figure 
5.  The  modes  or  shape  of  the  cylinders  circumferentially,  as  represented 
by  n,  are  believed  to  be  functions  of  the  shell  geometry,  material  type, 
and  the  applied  shock  loads.  If  analytical  methods  are  to  be  used  to 
predict  deformation  and  failure,  the  circumferential  mode  shape  must  be 
properly  determined  as  the  shell  stiffness  is  strongly  influenced  by  the 
number  of  buckles  generated. 


The  circumferential  mode  is  the  most  difficult  to  analyze  or  observe 
from  the  experiments.  Many  buckling  formulas  have  been  authored,  a major- 
ity of  which  were  determined  empirically.  Most  of  these  were  generated 
for  static  loads,  or  radical  loads  uniformly  applied  around  the  shell,  and 
were  not  applicable  for  the  conditions  of  transverse  blast  loads  as  applied 
in  these  exper iments.  rhe  experimental  data  gathered  here  tend  to  support 
theories  of  Greenspon  (.References  1 and  2)  which  suggest  that  a given 
cylindrical  shell  may  be  forced  into  a fundamental  buckling  or  collapse 
mode,  where  n might  be  considered  3 or  less,  or  a higher  buckling  mode  where 
n lias  some  large  value  that  is  a function  of  the  shell  length,  thickness, 
radius,  and  material  type.  Greenspon  further  suggests  the  use  of  Reynolds 
(Reference  ll)buckling  formula  to  determine  n.  Comparison  of  the  data  in 
Table  2 to  this  formulation  is  discussed  in  the  following  paragraphs,  along 
with  observed  trends  ;ind  deformation  hypotheses. 

The  data  points  in  Table  2 are  ordered  in  three  groups  of  constant 
L/D  values,  rhe  data  points  for  an  L/D  of  0.39  are  listed  for  decreasing 
a/h  values,  (Figures  b through  lb).  These  data  points  show  a decrease  in 
n as  the  thickness  increases.  This  trend  continues  for  L/D  values  of  0.89 
(Figures  17  through  28).  This  grouping,  however,  contains  three  data  points, 
19,  20,  and  22,  which  did  not  buckle,  but  appparently  deformed  in  the  collapse 
pattern.  The  data  for  L/D  values  of  1.89  show  five  collapse  patterns  at 
data  points  27,  28,  29,  31,  mid  33.  The  shells  of  this  group  that  buckled 
all  had  n values  between  10  mid  25.  This  indicates  that  for  L/D  values 
of  approximately  1 and  larger  the  higher  buckling  modes  are  less  dominant. 

The  shells  tested  in  this  group  are  shown  in  Figures  29  through  40.  The 
appearance  of  a collapse  mode  mid  buckling  mode  in  the  same  cylinder  size  is 
apparent  in  data  points  (18,  19),  (21,  22),  (30,  31),  and  (32,  33).  Bach 
set  of  data  .represents  cylindrical  shells  tested  with  the  same  L/D,  mid  a/h 
values  but  with  different  applied  loads.  This  suggests  that  for  a given 
shell  size  there  exists  a critical  load  that  determines  if  buckling  occurs 
or  a fundamental  collapse  pattern  is  formed.  Almroth  (Reference  12)  calcu- 
lated such  critical  loads  for  small  deflections  mid  assumed  elastic  defor- 
mation. Hxperimcntal  determination  of  dynamic  critical  load  values  would 
require  extensive  testing  and  it  appears  that  peak  pressure  or  impulse 
alone  is  not  a sufficient  description  of  the  load  to  use  as  critical  values 
to  determine  collapse  or  buckling.  To  vary  both  experimentally  appears 
impractical  outside  a shock  tunnel. 

The  data  points  in  Table  2,  which  indicated  buckling,  were  compared 
to  Reynold's  formula  Bquation  (4)  mid  compiled  in  Table  3. 

n«5ry[(l.b2b)(L/D)(a/h)ts  " A (4) 

Reynold's  approximation  for  n has  been  rewritten  in  terms  of  L/D  and  a/h 
for  ease  of  use,  mul  as  seen  in  Table  3 predicts  circumferential  mode  nvmbers 
reasonably  close  for  L/D  values  less  than  1.  As  the  L/D  value  goes  to  1.89 
the  prediction  formula  overpredicts  by  a factor  of  2.  The  use  of  such 
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formulas  is  limited  until  a satisfactory  model  for  separating  buckling  from 
collapse  patterns  is  derived.  Since  the  problem  appears  to  be  one  of  stabi- 
lity, the  determination  of  a critical  buckling  load  appears  the  most  att^' 
tive.  If  the  appearance  of  a buckling  pattern  can  be  established,  formulas 
independent  of  load,  such  as  Reynold's  formula  may  be  accurate  enough  for 
analysis  purposes. 

The  failure  of  the  shells  as  defined  in  this  report  was  considered  to 
be  material  rupture.  Briefly  discussed  in  the  preceding  paragraphs,  this 
was  observed  to  initiate  along  the  clamped  edge  at  the  0 - 0 point.  Column 
(8)  of  Table  2 indicates  fracture  or  nonfracture  of  the  shell  material. 

The  approximate  centerpoint  deflection  is  shown  in  column  (7).  As  would 
be  expected,  the  longer  shells  deflected  more  before  failure,  but  due  to 
lack  of  control  over  the  degree  of  fracture,  little  else  may  be  determined 
from  the  data.  It  is  believed,  however,  that  failure  could  be  related  to 
the  centerpoint  deflection,  and  for  a given  shell  there  exists  a critical 
centerpoint  deflection  for  failure. 


SECTION  V 


CONCLUSIONS 


Blast  loaded  cylinders,  deformed  plastically,  show  two  types  of  perma- 
nently deformed  patterns;  i.e.,  a collapse  pattern  and  a buckled  pattern. 
The  pattern  ultimately  formed  is  dependent  on  the  load  as  well  as  the  geo- 
metric and  material  properties  of  the  shell.  If  buckling  patterns  are 
formed,  the  circumferential  mode  number  established  as  permanent  set  in  the 
shell  appears  relatively  insensitive  to  the  load.  The  deformation  at  or 
near  failure  encompasses  only  25  to  30  percent  of  the  shell  circumference. 
Failure,  defined  here  as  fracture,  occurred  at  the  fixed  boundaries,  and 
appears  to  be  a simple  tensile  failure  as  a result  of  axial  strain. 
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Figure  5.  Typical  Shell  Failure 


Figure  31.  Cylinder  for  Data  Point  2e>,  Table  I! 
a/h  = 188  L/D  = 1.89 
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Figure  33.  Cylinder  tor  Data  Point  ~S,  I able  TI 


Figure  34.  Cylinder  for  Data  Point  29,  Table  II 
a/h  - 117  L/D  - 1.89 
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